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Background: Recent experiments on beta-delayed fission in the mercury-lead region and the discovery of asym- 
metric fission in 180 Hg 1 have stimulated theoretical interest in the mechanism of fission in heavy nuclei. 

Purpose: We study fission modes and fusion valleys in 180 Hg and 198 Hg to reveal the role of shell effects in 
pre-scission region and explain the experimentally observed fragment mass asymmetry and its variation with A. 

Methods: We use the self- consistent nuclear density functional theory employing Skyrme and Gogny energy 
density functionals. 

Results: The potential energy surfaces in multi-dimensional space of collective coordinates, including elongation, 
triaxiality, reflect ion- asymmetry, and necking, are calculated for 180 Hg and 198 Hg. The asymmetric fission valleys 
- well separated from fusion valleys associated with nearly spherical fragments - are found in in both cases. The 
density distributions at scission configurations are studied and related to the experimentally observed mass splits. 

Conclusions: The energy density functionals SkM* and D1S give a very consistent description of the fission 
process in 180 Hg and 198 Hg. We predict a transition from asymmetric fission in 180 Hg towards more symmetric 
distribution of fission fragments in 198 Hg. For 180 Hg, both models yield 100 Ru/ 80 Kr as the most probable split. 
For 198 Hg, the most likely split is 108 Ru/ 90 Kr in HFB-D1S and 110 Ru/ 88 Kr in HFB-SkM*. 

PACS numbers: 21.60.Jz, 24.75. +i, 27.70. +q 



I. INTRODUCTION 

The fission phenomenon is a magnificent example of a 
quantal large-amplitude collective motion during which 
the nucleus evolves in a multidimensional space repre- 
senting shapes with different geometries, often tunneling 
through a classically-forbidden region [2] . Understanding 
the fission process is crucial for many areas of science and 
technology. For instance, fission governs the existence 
of many transuranium elements, including the predicted 
long-lived super- heavy species. In nuclear astrophysics, 
fission influences the formation of heavy elements in a 
very high neutron density environment. Fission applica- 
tions are numerous. For instance, improved understand- 
ing of the fission process will enable scientists to enhance 
the safety and reliability of nuclear reactors. While in the 
past the design, construction, and operation of reactors 
were supported through empirical trials, the new phase in 
nuclear energy production is expected to rely heavily on 
advanced modeling and simulation capabilities utilizing 
massively parallel leadership-class computers. 

A comprehensive explanation of nuclear fission rooted 
in interactions between nucleons still eludes us, although 
self-consistent approaches based on the nuclear density 
functional theory (DFT) have recently demonstrated 
that a microscopic description has a potential for both 
qualitative and quantitative description of fission data [3]- 
[8]. A starting point in the adiabatic approach to fission 
is the capability to compute accurate multidimensional 
potential energy surfaces (PES), and use them to predict 
observables such as fission half-lives and fragment mass 
distributions. 



This work has been stimulated by recent experiments 
on beta-delayed fission in the mercury-lead region pQ and 
the discovery of asymmetric fission of the nucleus 180 Hg. 
Such an outcome has not been initially anticipated, as 
the symmetric fission channel involving two semi-magic 
90 Zr fragments was believed to dominate the process. It 
has been generally expected that the asymmetric fission 
is not important below 227 Th [9 , in particular in pre- 
actinide nuclei with high-lying saddle-point configura- 
tions that depend weakly on shell effects [TO] , Moreover, 
the data on mass distributions of fragments in the low- 
energy fission of nuclei with 187 < A < 213 have demon- 
strated the strong presence of the symmetric fission mode 
[TT] . In particular, the nucleus 198 Hg has been observed 
to exhibit a fairly broad mass distribution [TTJ [12] . 

The explanation of the asymmetric fission around 
180 Hg has been offered by the macroscopic-microscopic 
model [TJ [13] and its extension [14] in terms of an asym- 
metric fission pathway that is separated by a potential- 
energy ridge from the symmetric 90 Zr+ 90 Zr fusion valley. 
These results have emphasized the importance of shell ef- 
fects between fission saddle and scission in pre-actinide 
nuclei (see also Ref. [32]). In this work, we extend the 
theoretical analysis of Ref. [T] using the self-consistent 
nuclear DFT. We explain the transition from the asym- 
metric fission in 180 Hg [1 to a more symmetric situation 
in 198 Hg [11 in terms of shell effects. We compare the 
fission pathways in both nuclei and discuss the interplay 
between fission and fusion valleys. Finally, by studying 
density distributions of fragments, we demonstrate that 
scission configurations of 180 Hg and 198 Hg can be under- 
stood in terms of molecular structures. 
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FIG. 1. (Color online) Fission pathways for 180 Hg (top) and 198 Hg (bottom) as functions of the driving quadrupole moment, 
Q20 calculated with HFB-SkM* (left) and HFB-D1S (right). Competing fission and fusion valleys are indicated, together with 
the associate shapes. Since the sEF configuration sometimes corresponds to a ridge (rather than a valley), the corresponding 
curves do not continue in such cases. See text for details. 



II. MODEL 

The Skyrme-HFB calculations were carried out using 
the framework previously discussed in Refs. [5j [16j [17] 
based on the symmetry unrestricted DFT solver HFODD 
[T8] capable of breaking all self-consistent symmetries of 
nuclear mean fields on the way to fission. To solve a 
constrained nonlinear HFB problem precisely, we em- 
ployed the Augmented Lagrangian method [16 j. The nu- 
clear energy density functional was approximated by the 
SkM* functional [19] in the particle-hole channel and the 
density-dependent mixed pairing interaction [20] in the 
particle-particle channel. To truncate the quasi-particle 
space of HFB, we adopted the quasiparticle-cut-off value 
of 60MeV in the equivalent energy spectrum [21 . The 
pairing strengths were adjusted to reproduce the neutron 
and proton pairing gaps in 252 Fm [5]; the resulting values 
are V n0 = -268.9 MeVfm 3 and V p0 = -332.5 MeVfm 3 . 
The stretched harmonic oscillator (HO) basis of HFODD 
was composed of states having not more than TVo = 26 
quanta in either of the Cartesian directions, and not more 
than 1140 states in total. 

Our Gogny calculations follow the framework de- 
scribed in Refs. [3j [22] based on the axial Gogny-HFB 
solver [23] and D1S parameter set [24]. We used the 



stretched HO basis with N z = 22 HO shells along the 
symmetry axis and N± = 15 shells in the perpendicular 
direction. The oscillator length was adjusted at every 
calculation point. 

To find the optimum trajectories in a multidimensional 
collective space, we constrain the nuclear collective co- 
ordinates associated with the multipole moments Qa/x, 
by which we explore the main degrees of freedom re- 
lated to elongation (A/i = 20) and reflection- asymmetry 
(A/i = 30). In our symmetry unrestricted Skyrme- 
HFB calculations we also explore the effects of triaxiality 
(A/i = 22) and necking (A/i = 40). In our axial Gogny 
calculations, the scission configurations were studied by 
means of the neck coordinate Qn defined through the 
gaussian-type operator Qn = exp[— (z — zo) 2 /a 2 ] with 
a = 0.1 fm and zq chosen to describe the neck region 
(e.g., zq = 0.5 fm in 180 Hg). Qn describes the number 
of nucleons in a thin layer of thickness a perpendicular 
to the symmetry axis placed at z = Zq. Large values 
of Qn describe shapes with a thick neck. By decreas- 
ing Qn one can approach the scission line. To obtain a 
PES, constrained HFB equations are solved to minimize 
the total energy of the system at each point in the col- 
lective space. As demonstrated earlier [5] [6], exploring 
many collective coordinates makes it possible to identify 
saddle points [25 as the competing fission pathways are 
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usually well separated when studied in more than one 
dimension. 



III. COMPETING FISSION MODES 

Figure [l] shows the calculated fission pathways in 180 Hg 
and 198 Hg. Both models predict a fairly similar pat- 
tern. The reflection symmetric fission path associated 
with associated with elongated fragments (sEF) can be 
found for small deformations. The reflection- asymmetric 
path, corresponding to elongated fission fragments (aEF) 
of different masses and shapes, is branching away from 
the symmetric valley below Q20 = 100 b, and it passes 
through the mass asymmetric scission point (see Fig. [2] 
for a better view of aEF in the Q20 — Q30 plane). 
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FIG. 2. (Color online) PES for 180 Hg (top) and 198 Hg (bot- 
tom) in the plane of collective coordinates Q20 — Q30 in HFB- 
SkM* . The aEF fission pathway corresponding to asymmetric 
elongated fragments is marked. The difference between con- 
tour lines is 4MeV. The effects due to triaxiality, known to 
impact inner fission barriers in the actinides, are negligible 
here. 



Q20 ~ 200 b; they are associated with post-scission con- 
figurations in which the two fragments are well separated. 
Since the fission process is adiabatic, sFu and aFu are not 
expected to couple to aEF, which has a very distinct com- 
pound configuration exhibiting a pronounced neck even 
at Q20 ~ 250 b. 

The situation in 198 Hg is qualitatively similar but the 
mass asymmetry along aEF is significantly reduced, and 
the energy difference between aEF and sEF (correspond- 
ing to Q30 = in Fig. [2| is small, i.e., the PES is soft in 
the octupole direction as one approaches the saddle. In 
HFB-D1S calculations of Fig. [TJd), the energy balance 
between these two configurations is so fragile that a local 
transition from aEF to sEF, and back, is predicted at 
large deformations. This topography is consistent with a 
broad mass distribution of fission fragments observed in 
this nucleus. It is interesting to see that the magic struc- 
ture of 90 Zr manifests itself in a very low energy of sFu 
in 180 Hg ( 90 Zr+ 90 Zr) in both models. The asymmetric 
fusion valleys aFu become more favored in 198 Hg. 
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FIG. 3. (Color online) PES in HFB-D1S for 180 Hg (top) and 
198 Hg (bottom) in the (Q2CQ30) plane in the pre-scission 
region of aEF valley. The symmetric limit corresponds to 
Q30 = 0. The aEF valley is marked by a dashed line. Density 
profiles for various pre-scission configurations are indicated. 
The difference between contour lines is 0.5 MeV. Note differ- 
ent Q 3 o-scales in 180 Hg and 198 Hg plots. 



In 180 Hg, at large elongations, the aEF path is strongly 
favored over reflection-symmetric configurations (Q30 = 
0) associated with elongated fragments (sEF) and the 
symmetric compact fragment valley sCF - marked in 
Fig. [TJa). The fusion valleys, both symmetric (sFu) 
and asymmetric (aFu) appear very low in energy above 



IV. PRE-SCISSION CONFIGURATIONS 

The properties of fission fragments are governed by 
the nature of scission configurations at which a nucleus 
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FIG. 4. (Color online) (a) PES of 180 Hg in the (Q 20 ,Qiv) plane computed in HFB-D1S in the scission region of aEF. (b) 
Density distribution in 180 Hg close to scission at Q20 = 260 b, Qn — 0.3, and Q30 = 33.8 b 3 / 2 (marked by a star in panel (a)) 
compared to density distributions of 90 Zr (in its spherical ground state) and 72 Ge (in the excited deformed configuration with 
Q20 = 8 b). The density profiles for r± = (c), and for z = — 8fm and z — 10 fm (d) along the cuts marked by dotted lines in 
panel (b). 



splits [26 . The scission point is not precisely denned in 
the models yielding leptodermous densities. We assume 
that a scission configuration corresponds to a well-defined 
thin neck, and for greater elongations the neck decreases 
and the binding energy rapidly drops due to the Coulomb 
repulsion between the fragments. For the detailed anal- 
ysis, the hypersurface of scission points in the collective 
space has to be computed [7 j. To get more insight into 
the mass distributions of fissioning Hg nuclei, in Fig. [3] 
we show the topography of PES for 180 Hg and 198 Hg ob- 
tained in HFB-D1S for pre-scission configurations around 
the aEF valley. It is gratifying to see that the predictions 
of HFB-D1S and HFB-SkM* (shown in Fig. [2) are simi- 
lar. Namely, in both cases the aEF valley is separated by 
a ridge from the symmetric Q3o=0 line (sEF) for 180 Hg, 
while for 198 Hg the pre-scission surface is fairly soft in 
the Q30 direction, and the aEF pathway corresponds to 
much smaller mass asymmetries. 

The profiles of the density distribution for various con- 
figurations are also plotted in Fig. [3] The density profile 
corresponding to aEF (the most probable static scission 
point) in 180 Hg can be associated with a Ah / Al=99 j '81 
mass split, which is very consistent with the observed 
mass asymmetry of A#/A L =100(1)/80(1) (TJ. The min- 
imum of the neck joining the two pre- fragments is lo- 
cated at 0.6 fm from the center of mass. Density distri- 



butions for higher-lying pre-scission configurations shown 
in Fig. |3ja) are all similar, with the heavier fragment 
being nearly spherical and the lighter fragment elon- 
gated. When moving closer to the Q30 = line, we see 
that the neck is still pronounced and the fragments are 
elongated. That is, the symmetric configurations com- 
peting with aEF do not correspond to the two spheri- 
cal 90 Zr nuclei. A very similar situation is obtained in 
HFB-SkM*, where the predicted mass split at the static 
scission point is Ah /Al=101/79, and in the model of 
Ref. [13 : A H /A L « 103/77. 

The pre-scission shapes of 198 Hg are shown in Fig.j3^b) . 
The mass split along aEF is 108/90 and the minimum of 
the neck is located at about 1.7 fm from the center of 
mass. Interestingly, the heavier fragment is elongated 
while the lighter fragment is nearly spherical, i.e., this is 
exactly opposite to what has been predicted for 180 Hg. 
Again, our HFB-SkM* calculations yield a very consis- 
tent result, A^/A^=110/88, that is close to the mass 
split «1 11/87 of Ref. ^3]. 

Figure [4^a) shows the fission valley aEF for 180 Hg in 
HFB-D1S in the (Q20,Qn) plane as the scission point 
(small values of Qn) is gradually approached with in- 
creasing Q2o- For Q20 < 170b, the configuration sFu 
corresponding to two separated 90 Zr fragments lies above 
aEF (see also Fig. [ljb)). For Q20 > 170 b an energy 
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barrier appears between aEF and sFu. Moreover, as dis- 
cussed above, a transition from aEF to sFu is going to be 
strongly hindered by the very different intrinsic structure 
of these two configurations. 

we first considered the nuclei that (i) have the same 
N/Z ratio as the parent system; (ii) have mass numbers 
that reproduce the doubled mass of the outer part of the 
fragment situated outside the vertical line in Fig. [4^b); 
and (iii) have density distributions that match those 
of the fragments. The detailed analysis of the shape 
of 180 Hg in the near-scission configuration of aEF with 
Q20 = 260 b and Qn = 0.3 is shown in Fig.Qb-d). Since 
the nuclei 90 Zr and 72 Ge have the same N/Z ratio as 
180 Hg they are obvious candidates. The larger fragment 
of 180 Hg has a nearly spherical shape. This is consistent 
with the ground-state of semi-magic 90 Zr. The smaller 
fragment is strongly elongated. A deformed configura- 
tion of 72 Ge with Q2 = 8 b, at an excitation energy of 
3.4 MeV, fits the bill. By comparing the calculated den- 
sity profiles, we see that the pre-scission configuration of 
180 Hg can indeed be viewed as a molecular system con- 
sisting of spherical 90 Zr and deformed 72 Ge fragments 
connected by a thin neck with a slightly reduced den- 
sity that contains 8 protons and 10 neutrons. At scission 
the neck nucleons are shared between the two fragments. 
One likely split could be 100 Ru/ 80 Kr. 

Following similar analysis for 198 Hg presented in Fig.|5j 
we conclude that its near-scission aEF configuration can 
be viewed as a molecular system. It consists of 80 Ge in its 
nearly-spherical ground state with Q20 = 2.4 b, 98 Zr in 
a well-deformed (Q20 = 12 b) configuration, and a neck 
containing 8 protons and 12 neutrons. Therefore, one 
likely split could be 108 Ru/ 90 Kr. 

V. CONCLUSIONS 

Our self-consistent calculations based on the nuclear 
DFT with SkM* and D1S effective interactions give a 



very consistent description of the fission process in 180 Hg 
and 198 Hg. By considering several collective coordinates, 
we were able to follow static fission and fusion pathways 
in the configuration space. We confirm the findings of 
Ref. [1] that the asymmetric fission valley aEF is well sep- 
arated from fusion valleys associated with nearly spher- 
ical fragments. We conclude that the mass distribution 
of fission fragments in both nuclei is governed by shell 
structure of pre-scission configurations associated with 
molecular structures. In 180 Hg, both our models suggest 
100 Ru/ 80 Kr as the most probable split - a finding that is 
very consistent with experiment - and both predict sym- 
metric elongated configurations sEF to lie rather high 
in energy. The most likely split predicted for 198 Hg is 
108 Ru/ 90 Kr in HFB-D1S and 110 Ru/ 88 Kr in HFB-SkM*. 
Both models yield PES for this nucleus to be fairly soft in 
the Q30 direction in a pre-scission region, and this is ex- 
pected to result in an increased yield of nearly-symmetric 
partitions, and yield a very shallow, or even two-humped, 
structure seen experimentally [TTJ [12] . 
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FIG. 5. (Color online) Similar as in Fig.|Z]but for 198 Hg. Density distributions of 198 Hg in (b)-(d) correspond to Q20 = 330b, 
Qn = 0.33, and Q30 = 12.1b 3 / 2 . They are compared to those of 98 Zr (in the excited deformed configuration with Q20 = 12b) 
and 80 Ge (in its prolate ground state with Q20 = 2.4b). 
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